Inhomogeneous Fulde-Ferrell superfluidity in spin-orbit coupled atomic Fermi gases 
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Inhomogeneous superfluidity lies at the heart of many intriguing phenomena in quantum physics. 
It is believed to play a central role in unconventional organic or heavy-fermion superconductors, 
chiral quark matter, and neutron star glitches. However, so far even the simplest form of inhomo- 
geneous superfluidity, the Fulde-Ferrell (FF) pairing state with a single centre-of-mass momentum, 
is not conclusively observed due to the intrinsic complexibility of any realistic Fermi systems in 
nature. Here we theoretically predict that the controlled setting of ultracold fermionic atoms with 
synthetic spin-orbit coupling induced by a two-photon Raman process, demonstrated recently in 
cold-atom laboratories, provides a promising route to realize the long-sought FF superfluidity. At 
experimentally accessible low temperatures (i.e., 0.05Tf, where Tf is the Fermi temperature), the 
FF superfluid state dominates the phase diagram, in sharp contrast to the conventional case without 
spin-orbit coupling. We show that the finite centre-of-mass momentum carried by Cooper pairs is 
directly measurable via momentum-resolved radio-frequency spectroscopy. Our work opens the way 
to direct observation and characterization of inhomogeneous superfluidity. 

PACS numbers: 05.30.Fk, 03.75.Hh, 03.75.Ss, 67.85.-d 
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Interacting Fermi systems with imbalanced spin pop- 
ulations are quite ubiquitous in nature, with manifes- 
tations ranging from solid-state superconductors to as- 
trophysical objects Q. The spin-imbalance disrupts the 
Bardeen-Cooper-Schrieffer (BCS) mechanism of super- 
conductivity, where fermions of opposite spin and mo- 
mentum form Cooper pairs. As a result, an exotic super- 
conducting state, characterized by Cooper pairs with a fi- 
nite centre-of-mass momentum and spatially non-uniform 
order parameter, may occur, as predicted by Fulde and 
Ferrell nearly fifty years ago [2fl . More complicated inho- 
mogeneous pairing states are also possible with the in- 
clusion of more and more momenta for the spatial struc- 
ture of order parameter, following the idea by Larkin and 
Ovchinnikov (LO) [3|. These forms of inhomogeneous 
superfluidity, now referred to as FFLO states, have at- 
tracted tremendous theoretical and experimental efforts 
over the past five decades 0|. Remarkably, to date there 
is still no conclusive experimental evidence for FFLO su- 
perfluidity. In solid state systems such as the organic su- 
perconductor A-(BETS)2FeCi4 [j| and the heavy fermion 
superconductor CeCoIns [|| 0], the experimental diffi- 
culty arises from unavoidable disorder effects and or- 
bit/paramagnetic depairings close to the upper critical 
Zeeman field. 

An ultracold atomic Fermi gas has proven to be an 
ideal tabletop system for the pursuit of FFLO superflu- 
idity [1| . Although it is largely analogous to an electronic 
superconductor, the high controllability in interactions, 
spin-populations and purity leads to a number of unique 
experimental advances J8| . Indeed, following theoreti- 
cal proposals by Orso [9( and the present authors 

[13- 

12fl ) strong experimental evidence for the FFLO pairing 
has been observed in a one-dimensional spin-imbalanced 
Fermi gas of 6 Li atoms [l^|. In three dimensions (3D), 




FIG. 1: (color online) Phase diagram of a 3D spin-orbit cou- 
pled atomic Fermi gas at a broad Feshbach resonance and 
at the lowest experimentally attainable temperature 0.05Tf. 
The synthetic spin-orbit coupling is induced by two counter- 
propagating Raman laser beams with strength Q.r and detun- 
ing 8. Here we take the recoil momentum kn = Uf- A narrow 
LO wedge near £Ir ~ and 8 ~ 1.2Ef is not shown. By 
increasing 8, the Fermi cloud changes from a FF superfluid to 
a normal gas, via first-order (dashed line) and second-order 
(solid line) transitions at low and high respectively. The 
FF superfluid can be either gapped or gapless. These two 
phases are separated by the dot-dashed line. The color shows 
the magnitude of the centre-of-mass momentum of Cooper 
pairs for the FF superfluid, q/ftF- The BCS superfluid occurs 
at Qr — or 8 = only. 



unfortunately, the FFLO phase occupies only an ex- 
tremely small volume in parameter space [13. [TEl and 



I1LL2I. 



thus is impossible to observe experimentally 

In this Letter, we predict that the FF superfluid could 
be easily observed in a resonantly interacting 3D atomic 
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Fermi gas with synthetic spin-orbit coupling. This sys- 
tem was recently realized experimentally at Shanxi Uni- 
versity [lH and at Massachusetts Institute of Technol- 
ogy [ISfl) using a two-photon Raman process. Our main 
result is summarized in Fig. [TJ At experimentally ac- 
cessible temperatures (T — 0.05I» [lj] the FF state 
occupies a major part of the phase diagram. By tun- 
ing the strength (£Ir) and detuning (S) of two Raman 
laser beams, the centre-of-mass momentum of Cooper 
pairs q can be as large as the Fermi momentum kp. We 
propose that such a large FF momentum can be easily 
measured by momentum-resolved radio-frequency spec- 
troscopy (see Fig. HJ [2(| • We note that the emergence 
of a FF superfluid has also been predicted in a 3D atomic 
Fermi gas in the presence of 2D Rashba spin-orbit cou- 
pling (2l| or 3D isotropic spin-orbit coupling [22], or in 
a two-dimensional atomic Fermi gas with a two-photon 
Raman process [23]. These systems are yet to be realized 
experimentally. We also note that the possibility of the 
FF superfluidity in our setting has been discussed very 
recently by Vijay Shenoy 24 ]. 

We start by considering a 3D spin-orbit coupled spin- 
1/2 Fermi gas of 6 Li or 40 K atoms near a broad Feshbach 



resonance 
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[19]. Experimentally, the synthetic spin- 
orbit coupling is induced using two counter-propagating 
Raman laser beams (i.e., along the z-direction) that cou- 
ple the two different spin states, following the same 
scenario as in the NIST experiment for a 87 Rb Bose- 
Einstein condensate (BEC) [25|. The Raman process can 
be described by, (fl R /2) f dx[^\ (x) e l2kRZ ^! i (x) +H.c], 
where (x) is the creation field operator for atoms 
in the spin-state a, fl R is the coupling strength of Ra- 
man beams, and k R = 2ir/\ R is the recoil momentum 
determined by the wave length X R of the two beams. 
Thus, during the two-photon Raman process, a momen- 
tum of 2hk R is imparted to an atom while its spin is 
changed from |!) to |t). This creates a correlation be- 
tween the spin and orbital motion, which can be seen 
more clearly by taking a gauge transformation, ^(x) = 
e ifeflZ V>t( x ) and *j.(x) = e~ lkRZ iJj^(x.). Near a Feshbach 
resonance, the system may therefore be described by a 
single-channel model Hamiltonian H = J dx["Ho +Hint], 
where the single-particle part 



Ho — 



a + Afc, + 5/2 fl R /2 

& - \h - 8/2 











(1) 



and Hint = Uoip^ip^ipiipt ( x ) is the interaction Hamilto- 
nian that describes the contact interaction between two 
spin states with strength Uq. The interaction strength 
should be expressed in terms of the s-wave scattering 
length a B , i.e., 1/U a = m/(^irh 2 a s ) - l/V^rn/^k 2 ), 
which can be tuned precisely by sweeping an external 
magnetic field around the Feshbach resonance Q. Here 
V is the volume of the system. In the single-particle 
Hamiltonian (JTJ) , k z = —id z , £k = — h 2 V 2 /(2m) — fi af- 



ter dropping a constant recoil energy, and 6 is the two- 
photon detuning from the Raman resonance. For con- 
venience, we have defined a spin-orbit coupling constant 
A = H 2 k R /m. In the Shanxi experiment with 40 K atoms 
flSjj . the Fermi wavelength is about kp ~ 1.6/cr and 
the coupling strength Q, R ~ 1-5.Br ~ 0.6Ep, in units 
of the recoil energy E R = Y?k 2 R J [2m) or Fermi energy 
E F = h 2 k 2 F /{2m). 

For the model Hamiltonian Jl} , most previous theoret- 
ical studies focused on the case with (5 = and considered 
the standard BCS superfluidity [26], 27 1. In this work, we 



are interested in the FF pairing state with a single valued 
center-of-mass momentum. Our investigation is motived 
by the interesting finding that the two-body bound state 
at 5 7^ always acquires a finite center-of-mass momen- 
tum along the z-axis (28|. At the many-body level, we 
therefore anticipate that at finite detunings a FF super- 
fluid could be more favorable than a BCS superfluid [24(. 

By assuming a FF-like order parameter A(x) = 
— C/o (V^M^tM) = Ae lqz , we consider the mean- 
field decoupling of the interaction Hamiltonian, Hint — 
-[A(x)V>j(x)i/>J(x) + He] - A 2 /U . Then, using a 

Nambu spinor <I>(x) = [ijjf, ipi, tp^., ipJ] T , the total Hamil- 
tonian can be written in a compact form, % = 
(1/2) / dx$t( x )H BdG $( x )-l/A 2 /C/ + E k ek, where the 
Bogoliubov Hamiltonian 



HbcIG 



S£ n R /2 -A(x) 

0^/2 S~ A(x) 

A*(x) -S\ -Or/2 

-A*(x) -Qr/2 -SZ k 



(2) 



and iS^ = £k ± A/ ; z ± 5/2. It is straightforward to 
diagonalize the Bogoliubov Hamiltonian 'H.BdG^krii'^-) = 
Ek v &k v (x.) with quasiparticle wave-function $k^(x) = 
l/Vy e ik-x[ Uk ^ e+ i 9 V2 j u^e+^Z 2 , v^te-^/ 2 , w k „ie-^/2]T 
and quasiparticle energy E krj (rj = 1,2,3,4). The mean- 
field thermodynamic potential Q at temperature T is 
then given by 



Q 
V 



1 

2V 



E (^k+q/2 + Ck-q/2) ^ ^2 
k kr; 



k B T 
V 



EK 1 

k?7 



e -W]_i. (3) 

Uq 



Note that the summation over the quasiparticle energy 
must be restricted to E^ > because of an inherent 
particle-hole symmetry in the Nambu spinor representa- 
tion [29j. For a given set of parameters (i.e, the tempera- 
ture T, interaction strength \/kpa s etc.), different mean- 
field phases can be determined using the self-consistent 
stationary conditions: dfl/dA = 0, dil/dq = 0, as well 
as the conservation of total atom number, N = —dfl/dfi. 
At finite temperatures, the ground state has the lowest 
free energy F = Q + fiN. 
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FIG. 2: (color online) (a) Landscape of the thermodynamic 
potential, Cl(A,q) - 0(0,0), at tl R = 2E F and 8 = 0.6&Ef- 
The chemical potential is fixed to /i = — 0.471-Ef- The com- 
peting ground states include (i) a normal Fermi gas with 
A = 0; (ii) a fully paired BCS superfluid with A / and 
q = 0; and (iii) a finite momentum paired FF superfluid with 
A and q ^ 0. (b) The free energy of different competing 
states as a function of the detuning at Q.r = 2Ef- The inset 
shows the detuning dependence of the order parameter and 
momentum of the FF superfluid state. 



Without loss of generality, hereafter we consider 
the resonance case with a divergent scattering length 
l/k,Fa s = and set T = 0.05Tp, which is the lowest 
temperature that has been reported for Fermi atoms fl7l | . 
According to the typical number of atoms in experiments 
3 3|; we shall take kn = kp. corresponding to a di- 
mensionless spin-orbit coupling constant Xkp/Ep = 2. 
The Raman coupling strength f2# and detuning S can be 
easily tuned experimentally. Thus, we focus on the phase 
diagram as functions of £Ir and 6. 

In general, for any set of parameters there are three 
competing ground states that are stable against phase 
separation (i.e., d 2 fl/dA 2 > 0), as shown in Fig. HJa): 
normal gas (A = 0), BCS superfluid (A ^ and q = 0), 
and FF superfluid (A ^ and g ^ 0). Remarkably, 
in the presence of spin-orbit coupling (i.e., ^ 0) the 



FF superfluid is always more favorable in energy than 
the standard BCS pairing state at finite detunings (Fig. 
(Hb)). It is easy to check that the superfluid density of the 
BCS pairing state in the axial direction becomes negative 
(i.e., dft/dq < 0), signaling the instability towards a FF 
superfluid. Therefore, experimentally the Fermi cloud 
would always condense into a FF superfluid at finite Ra- 
man detunings. In Fig. [TJ we report a low-temperature 
phase diagram that could be directly observed in current 
experiments. The FF superfluid occupies the major part 
of the phase diagram. 

The greatly enlarged parameter space for a FF super- 
fluid can be qualitatively understood from the change of 
the Fermi surface due to spin-orbit coupling, as discussed 
in the pioneering work by Barzykin and Gor'kov (30| . 
In our case, the single-particle energy spectrum takes 
the form, E k ± = h 2 k 2 /2m± v /(fi H /2) 2 + (Xk z + 5/2) 2 , 
where "±" stand for two helicity bands. At large spin- 
orbit coupling (flu 3> Xkp, 6), we find 



£1 



h 2 (kl 



k± 



2m 



where q = 2kjid/flji. Thus, the Fermi surfaces remain 
approximately circular but the centers of the two Fermi 
surfaces are shifted by q/2 along the z-axis in opposite 
directions. In this situation, we may gain condensate 
energy by shifting the two Fermi surfaces by an amount 
±(g/2)e 2 , respectively 3l|. Here e z is the unit vector in 
z-axis. The fermionic pairing then occurs between single- 
particle states of k + (q/2)e z and k + (q/2)e z , leading 
to a FF order parameter that has a spatial distribution 
A(x) = Ae zqz . The direction of the FF momentum is 
uniquely fixed by the form of spin-orbit coupling (i.e., 
Xk z ) and its magnitude is roughly proportional to the 
detuning <5. 

At small spin-orbit coupling ~ 0, however, the 
above argument becomes invalid. Indeed, in the absence 
of spin-orbit coupling (f2ij = 0) the instability to inho- 
mogeneous superfluidity is driven solely by the detuning. 
The FF superfluid occurs in a narrow window (not shown 
in Fig. [J): S cs < S < 1.204E F , where S cs ~ 1.162^ 
is the Clogston limit for a unitary Fermi gas [15]. The 
direction of the FF momentum cannot be specified, in- 
dicating that a LO state with periodic stripe structure 
is more preferable [32]. In our calculations, we do not 
consider such a LO superfluid that exists in small wedge 
near ~ and 8 ~ 1.2Ep. Due to this simplification, 
at low Raman coupling strengths the FF superfluid does 
not transform continuously into a normal gas. 

The FF superfluid in Fig. [T] can have different 
topological structures in the quasiparticle excitation 
spectrum, as illustrated in Fig. [3l At small spin- 
orbit couplings (Fig. Eta)), the FF superfluid has a 
negligible center-of-mass momentum and therefore be- 
haves very similar to a standard gapped BCS super- 
fluid, while at large spin-orbit coupling (Fig. [Sib)), 



4 




(a) Q=0.6E F 



(b) Q=2.0E F 




FIG. 3: (color online) The excitation spectrum i?3(k) (in units 
of Ef) as functions of k± — ^Jk% + k'% and k z : (a) gapped 
spectrum at Qr = Q.QEf and (b) gapless spectrum at Q.r — 
2Ef- Here we take S = 0.5-Ef- Due to the particle-hole 
symmetry, E 2 (k) = -S 3 (-k) and J5i(k) = -S 4 (-k). The 
gapped branches -Ei(k) and i?4(k) are not shown. 



the sizable FF momentum features a gapless spectrum. 
The nodal points with zero-excitation energy are lo- 
cated at or close to the k z = axis, along which 
the dispersion relations take the form, E(k±,0) — 
±^/(Xq + S) 2 + ti 2 R /2 ± v/A 2 + {h 2 k\j2m + h 2 q 2 /8m) 2 . 
Apparently, these nodal points would appear when the ef- 
fective Zeeman field S 2 + Q 2 R /2 is sufficiently large. In 
Fig. [U the gapped and gapless FF superfluids are sep- 
arated by a dot-dashed line. The topological transition 
from gapped to gapless phases is continuous and could be 
revealed by thermodynamic measurements through the 
atomic compressibility, spin susceptibility and momen- 
tum distribution. 

We now consider how to experimentally detect a FF 
superfluid through momentum-resolved rf-spectroscopy 
[20l | . which is a cold-atom analog of the widely used 
angle-resolved photoemission spectroscopy in solid-state 
systems. In such a measurement, one uses a rf field with 
frequency lo to break a Cooper pair and transfer the spin- 
down atom to an un-occupied third hyperfine state |3). 
The momentum distribution of the transferred atoms is 
then measured by absorption imaging after a time-of- 
flight expansion. The rf Hamiltonian can be written as 
H rf (x j dyi[e~ ikRZ ipl (x) 4^ (x) +H.c], where t/>| (x) is 
the field operator which creates an atom in |3). The mo- 
mentum transfer kne z in T-L r f arises from the gauge trans- 
formation. For a weak rf drive, the number of transferred 
atoms can be calculated using linear response theory 



r(k,w) =A 



U 



(5) 

Here A u {k,uj) = J2 V l M k^|| 2 S(uj - E kv ) is the single- 
particle spectral function of spin-down atoms, £k = 
h 2 k 2 /(2m) - n, &ndf(x) = l/(e x / kBT + 1) is the Fermi- 
Dirac distribution function. For a FF superfluid, as 
Cooper pairs now carry a finite center-of-mass momen- 
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FIG. 4: (color online) Logarithmic contour plot of 
momentum-resolved rf spectroscopy: number of transferred 
atoms Y(k z ,u) at Q.r — 2Ef and at three detunings: (a) 
5 = and q = 0, (b) 8 = 0AE F and q ~ 0.1k F , and (c) 
5 = 0.8E and q ~ OMf. 



turn q = qe z , the transferred atoms acquire an overall 
momentum q/2. As a result, there is a q/2 shift in the 
transferred strength Eq. ([5]), which in principle could be 
experimentally measured. 

In Fig. 4, we report the momentum-resolved rf spec- 
troscopy along the z-direction T(k z , lu) = J2k ± ^O 4 , w) at 
£Ir = 2Ep on a logarithmic scale. Quite generally, there 
are two contributions to the spectroscopy, corresponding 
to two different final states: after a Cooper pair is bro- 
ken by a rf-photon, the remaining spin-up atoms can be 
in different helicity bands [33j|. These two contributions 
are well separated in the frequency domain, with peak 
positions indicated by the symbols "+" and "x", respec- 
tively. Interestingly, at finite detuning with a sizable FF 
momentum q, the peak positions of the two contributions 
are shifted roughly in opposite directions by an amount 
q/2. This provides clear evidence for observing the FF 
superfluid. 

In summary, we have shown that the condensate state 
of a spin-orbit coupled Fermi gas is the long-sought in- 
homogeneous Fulde-Ferrell superfluid. Our prediction 
can be readily examined in current experiments, where a 
three-dimensional spin-orbit coupled atomic Fermi gas is 
created using a two-photon Raman process [H,[l3. The 
finite center-of-mass momentum carried by the inhomo- 
geneous superfluid can be unambiguously measured by 
momentum-resolved radio-frequency spectroscopy. Our 
work complements relevant studies of solid-state systems, 
in which inhomogeneous superfluidity was predicted to 
be enhanced by Rashba spin-orbit interaction [30j, but 
was difficult to confirm experimentally. The controllable 
setting of a spin-orbit coupled atomic Fermi gas opens a 
new direction to explore the fascinating inhomogeneous 
superfluidity. 
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